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ABSTRACT 

We discuss the problem of observation of natural similarity in skeletal evolution of terres- 
trial mammals. Analysis is given by means of testing of the power scaling laws established 
in long bone allometry, which describe development of bones (of length L and diameter 
D) with body mass in terms of the growth exponents, e.g. A = d\ogL/d\ogD. The 
bone-size evolution scenario given three decades ago by McMahon was quiet explicit on the 
geometrical-shape and mechanical- force constraints that predicted A = 2/3. This remains 
too far from the mammalian allometric exponent A'^'^^p-' = 0.80 ± 0.2, recently revised by 
Christiansen, that is a chief puzzle in long bone allometry. We give therefore new insights 
into McMahon's constraints and report on the first observation of the critical-elastic-force, 
bending-deformation, muscle-induced mechanism that underlies the allometric law with es- 
timated A = 0.80 ± 0.3. This mechanism governs the bone-size evolution with avoiding 
skeletal fracture caused by muscle-induced peak stresses and is expected to be unique for 
small and large mammals. 

Keywords: allometric scaling , long bones, muscles, terrestrial mammals. 
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Observation of Mammalian Similarity through Allometric Scaling Laws 



by V.B. Kokshenev 



I. INTRODUCTION 

In general, biological laws do not follow from physical laws in a simple and direct way. 
Examples include Kleiber's allometric law known as the 3/4 power law that scales metabolic 
rates for animals and plants to their mass within the range of three order of magnitude. As 
shown by West et al. in Refs. |jl]J^ the observed metabolic rate scaling law arises from the 
interplay between geometric and physical constraints implicit, respectively, in space-filling 
fractal networks and energy dissipationB. Another famous 2/3 power law was proposed 
by McMahon 0] for scaling of longitudinal-to-transverse dimensions of animals and plants 
through physical description of geometric-shape and critical-force similarities noticed in 
their size evolution. Given in Ref. in explicit form, the geometrical- (cylindrical- volume)- 
shape and mechanical- (critical-elastic-buckling)-force constraints imposed on size evolution 
for animals and plants with their mass yielded the 2/3 power scaling law, along with the 1/4 
and 3/8 laws deduced respectively, for longitudinal and transverse linear dimensions. 
During almost three decades, McMahon's scaling laws have been a controversial subject of 
intensive study and debate. As a matter of fact, McMahon's description of the geometrical- 
shape and mechanical-force similarities was experimentally proved for terrestrial mammals 
neither in body allometry nor in long bone allometry [3,p| |l0tp!I] , |l2tp!3l , |l4tp^|l6| . Moreover, 
the most recent condemnation by Christiansen [0 states that no satisfactory explanation 



for any power-law scaling observed in mammalian allometry can be expected. 



*For simple explanations of Kleiber's allometric scaling, which is shown to originate from the 
general features of the networks irrespective of the geometrical and dynamical details, see Refs. 
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We will demonstrate that the failure of McMahon's constraints is due to the fact that the 
skeletal subsystem of animals is not mechanically isolated from their muscle subsystem, as 
was suggested in Ref. 0. Also, McMahon's hypothesis that the skeletal support of weight 
and fast locomotion of mammals is driven solely by a gravitation contradicts to up-to-date 
comprehension on a role of muscle fibers and tendons in formation of maximum skeletal 
stresses. We therefore revisit McMahon's evolution constraint equations in Sec. II within the 
context of their application to long bone allometry for terrestrial mammals. These equations 
are modified and generalized in view of the known experimental findings in muscle fiber 
allometry. Experimental testing of the two distinct critical-elastic-force mechanisms that 
govern evolution of mammalian bones is elaborated in Sec. III. Discussion and conclusions 
are given in Sec. IV. 



II. MCMAHON'S CONSTRAINTS IN LONG BONE ALLOMETRY 
A. Elastic Similarity Model Revisited 

Famous power laws by McMahon for scaling of linear dimensions of animals and 
plants was proposed within the framework of the so-called elastic similarity model (hereafter, 
ESM). Application of the ESM by McMahon to the case of mammalian bone allometry 
was based on the cylindric-shape correspondence that takes place between a given skeletal 
bone and a cylindrical beam. A bone sample was therefore geometrically approximated 
by a cylinder of diameter Dis and length Lis, where index i counts different bones and s 
indicates mammalian specie. The mechanical-force correspondence to the same rigid cylinder 
is justified by observation of the universal (specie-independent) bone-stress safety factors. 
These are given by ratio (about 3) of yield stress to peak stress. Exploration of such a kind of 
mechanical correspondence by McMahon gave rise to the maximum- (elastic-buckling)-force 
constraint imposed on volume-size evolution of a given bone. 

More specifically, the ESM is based on the fact that the mechanical failure of a bone 
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is prevented through its hnear dimensions Dig and Lj^, adjusted to bear critical buckling 
deformations, related to peak stresses through the maximum elastic forces: F^^^^^ = Pj^^"^]. 
The latter is describes an elastic instability caused by critical bending deformations that 
was specified by the Euler critical estimate -Ffc^cM = n'^EI/L'^ for a given cylinder (of length 
L and of diameter D, with the moment of inertia / = ttD^/GA and the elastic modulus E, 
see e.g. Cap. IV in Ref. ^^). Thus the ESM constraint equations attributed by McMahon 
to the cylindric-shape and elastic-force similar skeletal bones can be introduced through (a) 
the elastic-buckling critical force F^^'^'^*^ and (6) the cylindric-bone volume V^g"'^'^\ namely 

- (la) 

V,. = DlU, . (lb) 

In long-bone allometry, the observation of evolution of limb bones across mammalian 
species is discussed though the bone-size linear-dimension scaling to body mass Ms- This is 
given in terms of the bone-diameter and the bone-length allometric exponents, respectively, 
di and /j, or of the i-bone-dimension growth exponents, introduced by the following scaling 
differential relations, namely 

dXogDis dXogLis . . 

dlogM^ ctlogM^ 
^ ^ d\og Lis ^ h_ 
dlogDis ~ di 

The reduced dimension exponent Aj, related to the longitudinal-to-transverse scaling, is 
also defined. As seen, Eqs.(2),(3) are equivalent to the corresponding differential equa- 
tions dDis/dMg = diDis/Ms, etc., which solutions are commonly derived in bone allometry 
through regression equations Dis = CisM'^\ where M is treated as an external mammalian 
parameter and Cis are constants. A notable feature of the introduced scaling differential 
relations is independence of the i-bone exponents on mammalian specie s. This corrobo- 
rates the bone allometry observations and Eqs.(2),(3) are therefore treated as the allometric 
scaling laws. This implies a universal fashion in evolution of any linear dimension of bones, 
as well as bone volume Vis = Df^Lis, with body mass that in a certain way reflects similarity 




of mammals with their size evolution. With taking into account that pVis = Mis (p is bone 
density), and adopting additionally McMahon's hypotheses that (a) effective skeletal growth 
is driven by gravitation, i.e., F^^"^^ ~ gMis {g is the gravity constant) and that (6) bone 
mass Mis linearly scales to body mass Ms, the following i-bone-evolution equations, namely 



(4) 



result from, respectively, Eqs.(la) and (lb) with the help of the scaling relations given in 
Eq.(2). In turn, Eqs.(4) and (3) provide the well known ESM predictions: dg^"^'^'^ = 3/8, 
j{bucki) _ ]^^4^ g^j^^ y{hucki) _ 2/3^ including a trivial isometric solution do = = 1/3 and 
Ao = 1. As mentioned in the Introduction, these predictions were not experimentally proved 
even when a statistical dispersion of allometric data was taken into account (for recent 
criticism of the ESM predictions for the allometric exponents d, I and A see analyses given 
in Table 5 in Refs. and Table 3 in Ref. respectively). 



B. Elastic-Buckling- Force Criterium 

Skeletal evolution of animals cannot be studied independently of their muscle fibers and 
tendons. Moreover, the peak skeletal stresses are generated rather by muscle contractions 
than by gravitation. These both statements follow from studies of muscle design and bone 
strains during locomotion []TR|IU|,pU|,pT| . We infer therefore that the maximum elastic forces 



exerted by long bones are originated from the maximum muscle forces, i.e., P'^f^^f' = Fj^^^K 
The same studies provide strong evidence that the maximum muscle stresses are independent 
of body mass, and thus F^^^^A^^^] oc M°, where A^^^^] is the maximum cross-section 
area of muscle fibers. The critical- force constraint, justified through the aforementioned 
bone-stress safety factors, can be therefore formally introduced into consideration by the 
"overall-bone" critical-force exponent ac, namely 

_ dlogFj-j) _ _ d\ogA(Z7} 
'~ dhgMs ~ dlogM 



and the corresponding critical muscle-area exponent Ocm- These should be distinguished 
from the exponents 

C^log^if"*'' , d\0g Amuse /„v 

= ^logM. "™ = rflogM 

where F^^"^^ is given in Eq.(la). The muscle-area exponent is known in muscle allometry 
10| , |20| , ^ as the muscle-fiber, cross-section-area exponent and can be exemplified by data 
a^^P-* = 0.69 — 0.91 derived by Pollock and Shadwick for four distinct groups of muscles 
in mammalian hindlimbs (see Fig. 3 in Ref. ||22|). The maximum muscle force is commonly 
associated [l^ with the leg group of muscles of animals, i.e., A^^^^] = A^^^}^. With adopting 



of the latter in Eq.(5), the "leg-muscle" cntzca/ exponent a'f^^ = 0.81 — 0.83 was obtained 
on the bases of data [|1^ for six groups of mammalian leg muscles by Alexander et al.) 



and reported by Pollock and Shadwick in Ref. |22]. As seen, the means a^^^-* = 0.80 and 



^cm^^ = 0.82 are different but not distinguished within the experimental error. 

In order to establish the critical muscle-area exponent a^^^^ defined in Eq.(5), we have 
reanalyzed the experimental data by Pollock and Shadwick on muscle fiber area Amuse in 
hindlimbs of 35 quadrupedal mammalians as a function of body mass represented in Fig.l 
from Fig. 3 in Ref. [^]. In general, the gastrocnemius group of muscles (shown by diamonds 
in Fig.l for points adjusted |^2[ with a^-* = 0.77), unlike the common digital extensors 



group (shown by crosses for points with = 0.69), plays a principal role in 

formation of maximum bone stresses. This is due to the fact that for small and large species 
of animals > A'^\ Other two groups of muscles exhibit a crossover from the almost 
isotropic evolution with oq = 2/3 to somewhat given by the exponent acm and driven by the 
maximum muscle areas controlled by the gastrocnemius group, i.e. by Acm, established by 
maximum points of A^^^^^. First the highest 5 points (shown by arrows in Fig.l) are fitted 
by A^J^ = 428 * M°'^°. Next nearest-neighbors (16 points indicated as solid symbols below 
the line Acm in Fig.l) are fitted with Afi = 30A*M^-^\ Re gression elaborated within all the 
field of a maximum muscle area, defined by the highest 21points, provides = 0.82±0.01 
(with correlation coefficient r = 0.996). Remarkably, this finding matches well the aforegiven 



data for the " leg- muscle" exponent by Alexander et al. reported in Ref. ||22| and can be 
therefore treated clS cl reliable data. 

Eqs.(5) and (6) provide the following definition for the "overall-bone" averaged expo- 
nents, namely 

a, =< a,, >=-iz = with ai^^P) = 0.82 ± 0.01, (7) 

^ 1=1 

where summation is limited by bones which do play a principal role in effective support and 
fast locomotion of body mass of animals. Eq.(7) can be also treated as an extension of a 
similar definition of the mammalian principal-bone-averaged exponents d, I and A introduced 
with the help of Eqs.(2),(3), e.g., d =< di >. Thereby, revision of McMahon's a-hypothesis 
provides a new a-constraint equation imposed on the exponents: Ad — 21 = ac. 

In view of the fact that neither skeletal mass [|1^] nor bone mass |jT^ are linear with 
mammalian body mass, McMahon's revised 6-constraint equation given in Eq.(4) for i-bone 
is also modified as 2di + k = h^. Here the i-bone-mass exponent, namely 

dlogMis 

is introduced through the relevant power scaling law. Thus, McMahon's critical-force and 
cylindric-shape constraints given in Eq.(4) result in the following modified ESM constraints, 
namely 

4:d-2l = ttc, 

(9) 

2d + l = b. 

In turn, this yields new predictions for the mammalian overall-bone dimension and reduced- 
dimension exponents, or the elastic- buckling- criterium predictions, namely 

^ibucki) ^ ^c + 26 ^ ^(,,,,,) ^ 2b-a^ ^^^^ 
8 4 

^{bucki) ^ g ^ b_i ^ _2 _ 

aci + 2h 

The latter prediction follows from the definition for the reduced-dimension exponent A, = 
U/di given in Eq.(3) and presented here in the form Aj = bi/di — 2, with the help of the 
6-constraint equation. 



C. Elastic-Bending-Force Criterium 



After Alexander et al. it has been widely recognized (for recent references see Ref. 



|Tq| ) that the elastic bending deformations play a crucial role in the overall peak stresses 
of long bones instead of a simple axial compression discussed by McMahon in terms of 
the critical buckling deformations. The corresponding critical force Fj:^^"^^ = -Fften? applied 
normally to the bone before fracture was already discussed in long-bone allometry in Refs. 
p!^ , |20[] . In view of the elastic nature common for both kind of deformations, the force 
^bend ^ Certain way extends the ESM given in Eq.(l) for the case of the bending critical 
deformations, namely 

pDlLis = Mis. (12b) 

Straightforward application of the scaling differential relations introduced in Eqs.(2),(3), 
with accounting of the critical-force and the bone-mass growth exponents given in, respec- 
tively, Eqs.(5),(8) and (9), results in the following new constraint equations: 

(13) 

2d + l = b. 

This provides the elastic-bending criterium expressed in terms of the following predictions 

for the mammalian bone-dimension growth exponents, namely 

^ibend) ^ ^c + b ^ ^(,,„,) ^ 36-2a, 
5 5 

ttci + bi 

Notably that both the elastic-force criteria given in Eqs.(lO) and (14) are consistent with 
the isometric solution (do = = 1/3 and Aq = 1), which is found under conditions that the 
mammalian muscle-area subsystem develops isometrically (oq = 2/3) and independently of 
the skeletal subsystem [b^ = 1). The observed allometric scaling laws with d^'^^P) > 0.33, 
/(cxp) ^ 0.33, and X^"^^^^ < 1 corroborate that this simplified geometric scenario is avoided 
by the nature. 

8 



III. OBSERVATION OF BONE EVOLUTION SIMILARITIES THROUGH 
EXPERIMENTAL TESTING OF THE CONSTRAINT EQUATIONS 



All predictions given by the original ESM and the revised ESM are analyzed in the 
bone growth diagram in Fig. 2. As seen, the available experimental data matches neither the 
isometric nor the original ESM solutions (shown by crosses), even in case when dispersion 
effects of the experimental data (shown by error bars) are taken into account. Note that 
this large dispersion is not caused by error measurements of bone dimensions or body mass 
of animals, but is resulted from a large phylogenetic spectrum of terrestrial mammalsi. 
Unlike the case of the pioneer data ^ by Alexander et ai, all species which have multiple 
specimens, were additionally averaged [|l5l within a certain mammalian subfamily before 
to be documented. The most accurate allometric data with the systematically reduced 



phylogenetic statistical error was given ||14| , p!5| , p!6| by Christiansen. 

Predictions of the modified ESM and the extended ESM are shown in Fig. 2 by the shaded 
areas, which correspond, respectively, to Eqs.(9) and (13) estimated with account of the 
reliable domain for the critical-force exponent a^^^^^ = 0.81 — 0.83 and of that for the bone- 
mass exponent fo^'^^P) = 1.0 — 1.1 (that approximately covers error scatter of the experimental 
data on 6-'^^^'' (given in Table 2 in Ref. [0). The shaded areas indicate the critical- force 
constraints given by the a-constraint lines 4d — 21 = 0.82 and 3d — I = 0.82 extended by 
cylindric-volume constraints implicit in the form of the elastic-buckling-force and elastic- 
bending-force criteria, respectively. As seen from Fig.2, the elastic-buckling criterium seems 
to be observable within the range of the unreduced phylogenetic statistical error. After 
reduction of this error, only the elastic- bending criterium corroborates. 



Besides the case of the 6- long-bone-averaged allometric data |14] given in Fig.2 for the 



^In fact, there exist a certain error due to deviation of bone shape from the ideal cylinder. Also, 
not all body mass were really measured but taken as an average from the literature data (see 
discussion in Ref. [p^). 
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one-scale least-square regression {LSR), we have also elaborated analysis of the double set 
of the allometric exponents (taken from Table 5 in Ref. ||T^) derived within the two-scale 
regressions made for small (M < 50kg) and large (M > 50kg) mammals. But no definitive 
conclusions on domination of any elastic-force criteria can be inferred. Indeed, in the case 
of the overall- (6 -bone)-average analysis, the data for small and large animals is far to be 
fitted by any of the dashed areas in Fig. 2. When the ulna and the fibula are excluded, the 
principal- (4-bone)-averaged LSR data justifies the elastic-bending and the elastic-buckling 
criteria for the cases of small and large mammals, respectively. However, unlike the case of 
the one-scale data, experimental accuracy of the two-scale analysis is marginal that makes 
doubtful any inference on observation of both the distinct critical-force constraints. We 
have therefore restricted our analysis by one-scale allometric data for the four principal 
mammalian long bones listed in Table 1. 

First, we check a self- consistency of experimental data on the dimension ( df^^^ 
) [|I^ and reduced-dimension (A^'^^'' ) |[T^ allometric exponents obtained independently 
and presented in first and second columns of Table 1, respectively. As seen, when the 
bone-averaged data is compared between the two regression methods, it obeys the relation 

^(exp)/^(exp) ^ ^(oxp) ^j^j^ 

accuracy that is much higher than that for the case of partial i- 
bone relations df^^^ /I'f^^^ = Xf^^^ compared within the same method. Then, the geometrical 
mammalian similarity is tested on the basis of the 6-constraint equation 2d[^^'^^ -{-if'^'^^ = hf^^^ 
in second and third columns of Table 1. Again, the cylindric-shape constraints, given in terms 
of the bone-averaged data, are justified! with a good precision. We infer that observation of 
the mammalian similarity through the allometric power laws can be realized only "on the 
average" , but not for a given type of " mammalian" bone as it widely adopted in allometric 
studies. Examples are analyses of the original ESM predictions elaborated for a given z-bone. 



■^Exclusion should be given for the case of the exponent h* , which data obtained by the square 
regression (LSR) method is not available. 
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instead of the overall-bone data, and given in Table 5 in Ref. [|T4|, Table 3 in Ref. |T5|, and 



Table 3.11 in Ref. ^ 



The problem of validation of the bone-evolution a-constraint equation for the case of the 



bending loads, i.e., 3d — I = a, where a is treated as a free parameter, was first discussed ||20 
by Selker and Carter in terms of the bone strength index. On the basis of the mammalian 
data fill] by Biewener (shown in Fig. 2) and their own data for artiodactyls, the overall-bone- 



averaged equation Sd^'^^'^^ _/(exp) _ provided |^ estimates a = 0.77 and 0.82, respectively. 



By generalization of these findings to the overall mammalian case, allometric exponent 
a^^P-* = 0.77 — 0.82 was adopted pO| for testing of the bending or torsion deformations in 



mammalian long bones due to muscle contractions. The same analysis made on the basis 
of other available in the biological literature allometric data, including the particular case 
of birds!, has been recently given by Garcia (see Table 3.12 in Ref. |2^). As the result, 
the allometric muscle-area exponent was sugg ested a(^^P) = 0.77 - 0.83, with the mean 
a!^^^ = 0.80, as a suitable data for experimental testing of the bending-force constraint 
equation (see analysis in Table 3.11 in Ref. This suggestion is not true. 



Indeed, as follows from the pioneer work ^ by McMahon, revisited in the previous 
section, the force-constraint equation is driven by the critical force and therefore given as 
3d — I = Qc where the critical-force exponent, according to Eqs.(5),(7), is established by 
the data on maximum-muscle-area allometry, i.e., = a^cm'^ — O-^l ~ 0.83, with a^^^-* = 
0.82 that should be distinguished from the suggested ||2^ data on a^"^^ = 0.80. We have 



therefore reconsidered analysis given in Table 3.11 in Ref. ||2^ and found that no definite 
conclusions can be made on validation0 on the principal-bone averaged equation Sd^^^"^^ — 
^(cxp) _ q(^p) on the bases of the two-scale RMA and LSR data [|^. Conversely, the critical- 
bending-force constraint equation 3d — I = ac is strongly supported by the one-scale data 



^Application of the ESM for birds remains questionable. 

** Again, the marginal estimate acm = 0.829 is obtained in the case of the small-animal LSR data. 
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1^ by Christiansen deduced through both the different {LSR and RMA) regressions. This 
follows from the bone-dimension predictions = 0.82 and 0.83 (obtained I^S], respectively, 
for both the methods with the help of data given in first column in Table 1). 

In the current study we put emphasis on observation of the mammalian similarity through 
the critical muscle allometry exponent a^^m'^ established in Fig.l along with the one-scale 
long-bone allometric data on the reduced- dimension exponent A obtained in Ref. [[15 



Therefore, we have reformulated the elastic-buckling and the elastic-bending criteria given 
in Eqs.(ll) and (15) in terms of the observable Aj . This provides the following predictions 
for the critical-force exponents, namely 

(bucki) ^ 2 < ^^6, > and a^^^"'^) =< ^—^^k > (16) 

obtained with the help of Eq.(7) and estimated in last column of Fig.l. As seen, the elastic- 
force bone-^txcMmg'-deformation mechanism, proposed by McMahon in Ref. suggests an 
estimate a(^"^'=0 = o.90 for both the regression methods that is not justified by Eq.(7). 
In contrast, the elastic- force bone-^endrng^deformation mechanism predicted by af'^'^'^^ = 
0.81 (and 0.83) within the LSR (and RMA regression) methods is proved by the reduced- 
dimension and self-consistent linear- dimension allometry data reported by Christiansen in 



Refs. [[Tsl and Ref. |jT^, respectively. Again, analysis, similar to that given in Table 1, 



extended to the case of the two-scale principal-bone data 0,|l^], corroborates the same 
bending-force mechanism only for the data derived for small mammals within the LSR 
method. No conclusions can be inferred in the case of small mammals treated by the RMA 
regression, as well as in the case of large mammals treated by both the methods. 

IV. DISCUSSION AND CONCLUSIONS 

We have discussed the problem of observation of natural similarity in evolution of terres- 
trial mammals. Testing of the two conceivable underlying mechanisms that drive the bone 
size development with body mass of animals is given on the basis of experimental data on 
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the reduced dimension (Aj = U/ di), longitudinal dimension and transverse dimension 
(dj) allometric exponents established through the scaling laws commonly discussed in the 
mammalian long-bone allometry. 

Since Galilei it was repeatedly recognized that the isometric skeletal evolution prescribed 
by the overall-bone exponent Aq = 1 is not observed in the nature because the small mam- 
mals are not geometrically overbuilt and the large species do not operate very close to 
their mechanical failure limit, that it would be expected from the isometric scenario. This 
figurative, widely cited description given by Biewener [^] is in agreement with the simpli- 
fied version (Aq = 1) of the more sophisticated scenario (^x^^^^''^^ = 0.667) proposed by 
McMahon. Within the ESM, the mammalian similarity was introduced on the basis of 
realistic geometrical-shape and mechanical-force correspondence that takes place between a 
given skeletal bone and a rigid cylinder. As mentioned in the Introduction and illustrated in 
Fig. 2, predictions for evolution of bone dimensions with body mass given by the ESM were 
disapproved in long bone allometry. This, in particular, implies that the ESM exponent 
was not justified in long-bone allometric experiments, including the most 
systematic data with X^^^^^ = 0.78 — 0.82 (that follows from Table 1 as the mean between 
the LSR and RMA bone-averaged data). 

A good deal effort has been undertaking in long bone allometry to learn experimen- 
tal conditions for observation of the critical- force (elastic- ^wc/s/mg^deformation and gravity- 
induced) mechanism suggested by McMahon for explanation of anatomical adaptation of 
skeletal bones through their linear dimensions. The first objection by Economos was as 
follows. McMahon's mechanical-failure mechanism should not be expected as a unique for 
all species, but would more suitable for large mammals. This stimulated a careful search for 
additional scaling laws related to small and large mammals. Such kind of new scaling laws 
were established in terms of the double sets of allometric exponents introduced |T^,|T5|JT^ 
by Christiansen through the two-scale regressions distinguished by Mc = 50kg adopted as 
a boundary mass common for all species. Furthermore, it was speculated that the revealed 
inadequate description of the scaling laws is due to inaccuracy of the methods of regres- 
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sion and, as a result, the RMA regression was suggested as well-chosen instead of the 
traditional LSR. The second objection 0] by Economos refers to the linearity of the loga- 
rithmic scaling laws given in Eq.(2), which was not expected to be a sole across the three 
order of magnitudes of body mass. Experimental verification of this idea by Christiansen 



revealed WM that an application of the polynomial type of regressions in bone allometry 
does not improve the correlations established within the traditional linear logarithmic scal- 
ing. Finally, thorough numerical analyses ||l^ , p!5| of the reasons of the ESM failure brought 



Christiansen to a conclusion that "many factors contribute to maintaining skeletal stress 
at uniform level", including the factor of bending-deformation-induced stresses, which are 



more important [|I4| than the bone stresses illuminated ^ by McMahon. 

We have demonstrated how the factors of muscle fiber contractions, bone mass evolution, 
and of bending bone deformations can be incorporated into the ESM constraint equations. 
As a result, the modified (by bone- mass and muscle-contraction factors) ESM becomes ob- 
servable (see shaded area that extends a-buckling line in Fig. 2) under condition that the 
unreduced phylogenetic statistical dispersion of the allometric data P JTT|JTB[] is taken into ac- 
count. Otherwise, the extended (additionally by the 6en(img^deformation factor) ESM agrees 
with experiment (see Fig. 2). Another analysis (given in Table 1) yields the observation of 
the mammalian similarity within the principal-long-hone allometric data [ll4| , |T5i p!6| with sys- 



tematically reduced statistical error by Christiansen. As demonstrated, this observation is 
realized in terms of the bone- afera^'ed allometric exponents, restricted by the principal hones 
that are involved into the evolution-constraint equations. Example is the volume-constraint 
equation, which should be valid for any conceivable bone-evolution mechanism. As seen from 
analysis given in columns 2 and 3 in Table 1, the volume-constraint equation is much better 
"observed" in the "bone-averaged" form 2d^'^^'^^ + /('^^p) = 5('^^p) than in the form presented 
for a given z-bone. We guess that the observation of the geometric-shape similarity through 
experimental justification of exact Eq.(12b) should depend on neither the number of scales 
nor on methods chosen for regression of the bone-dimension allometric data. Our additional 
verification of the cylindric-shape similarity given on the basis of the two-scale principal- 
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long-bone allometric data (taken from Table 3 in Ref. [|1^]) corroborates this statement for 
both small and large mammals. We infer therefore that both the methods and both the 
scales are equivalent in observation of the "bone-averaged" geometrical-shape mammalian 
similarity, at least for the principalH bones. 

As highlighted by Christiansen, the principal bones play a crucial role in primarily sup- 



port the body mass. He noted |15[ that greatly reduced ulna and too thin fibula do not play 
of much importance in support of body mass. They are therefore not suitable for testing 
of the critical-force constraints and should be excluded from the principal bone set. As 
follows from our analyses given in Fig.l, qualitatively the same should be referred to some 
muscle fiber groups such as common digital extensors which eventually do not produce peak 



bone stresses. Meanwhile, as seen from Fig. 9 in Ref. [0] and Fig.l, tibia and plantaris 



show a crossover behavior between principal and non-principal sets, where the principle 
bone-set and the principle muscle-group set are presented by small and large mammals, 
respectively. As the reliable cntzca/ principal- muscle-area exponent (= etc), which enters 
the critical-force a-equations, the data a^^m'^ — 0-^2 ± 0.01 is proposed in Eq.(7). This find- 
ing is derived in Fig.l from the maximum muscle-area gastrocnemius group of leg muscles 
and should be distinguished from the overall muscle-area data a(^^P) = 0.80 ± 0.03 that was 
groundlessly used, instead of a^^P\ in establishing of experimental validation [^,^ of the 



critical- ^endzng^force constraint 3(i*^'^^P-' — /('^^p) = a^^P-*. As shown, this equation, unlike the 
case of the ciitical-buckling-foice constraint Ad^"^^^^ — 2/(°^p) = a^^P-* related to the original 
ESM, is observable directly and indirectly through, respectively, the a-constraint equation 
and Eq.(16) (analyzed in the last column of Table 1). Again, we infer that the observation 
of the bending-ioYce criterium does not depend on the method chosen within the one-scale 
regression. 



^I^Extended statistical analysis of both the constraint equations, with including all available bone 



allometric data will be discussed elsewhere |25| 
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This is not the case for the two-scale data on bone-dimension allometric exponents re- 
ported ||T^ by Christiansen. Indeed, as follows from our many-sided analysis, the elastic- 
bending criterium is definitely supported for small and large mammals within the (principal- 
bone-averaged) LSR data and RMA, respectively. With accounting of the observation of the 
same criterium though the one-scale (principal-bone-averaged) LSR data, we see that corre- 



lations established by the traditional LSR method, unlike suggestion given in Ref. 0], show 
their self-consistency. But no certain conclusions can be inferred within the observation win- 
dows for small and large mammals in the cases of, respectively, RMA regression and LSR. 
We guess that the revealed discrepancy of the two equal in rights regression methods signals 
on failure of definition of the observation windows employed for the analysis of the critical- 
force constraints. In other words, unlike the case of the cylindric-shape similarity, these 
windows are not expected to be universal for observation of the elastic-force mechanisms, 
and cannot be therefore introduced by the unique boundary mass Mc- 

Thereby we have demonstrated that the mammalian similarity, observable through ex- 
perimental validation of the bone-evolution constraint equations, is described in terms of the 
one-scale principal-hone- averaged characteristics, which show independence on the regres- 
sion methods. Within this context, the observed in nature long-bone mammalian evolution 
can be described through longitudinal-to-transverse bone-dimension scaling law, with the 
aforegiven "method-averaged" exponent X^^^p^ = 0.80 ±0.02. Assuming a high enough accu- 
racy for the i-bone experimental data on the exponents a^^i^"^^ ^'f'^^\ both the discussed 
evolution-mechanism criteria are approximated in the following forms, namely 

yibuckl) ^ ^ 2& — acm y(bend) ^ 3& — 2acm 

26 + acm h + acm 

that follows from Eqs.(ll) and (15), respectively. With accounting of a'f^^ = 0.81 — 0.83, 
and adopting for the bone-mass mammalian allometric exponent ^(''^p) = 1.03 — 1.06 (see 
column 3 in Table 1) one has the following reduced-error estimates for the longitudinal-to- 
transverse scaling allometric exponents: 

^{buckl) ^ Q g7 ^ Q Q2 )^{hend) ^ Q gQ ^ Q ^j^j^ ^(exp) _ q gg ^ Q Q2. (18) 
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One can see that solely the elastic-bending criterium is validated. This implies corroboration 
the bone evolution mechanism, which provides avoidance of mechanical failure of mammalian 
bones caused by critical elastic bending deformations induced by maximum- area muscle 
contractions achieved in long bones during peak stresses. 

Prom the physical point of view, the fact that the bending (but not buckling) elastic 
deformations are crucial for mechanical failure of long rigid bones should be expected, under 
condition that the inequality Lis S> Dig (but not Lis ~ Dis) is fulfilled for animals of arbi- 
trary mass. Meanwhile, this fact was not corroborated in the one-scale long-bone allometry, 
and we therefore report on the first observation of the bending-critical-force bone-evolution 
mechanism, which is suggested to be universal regardless of small and large mammals. Fi- 
nally, after McMahon, we have demonstrated how the scaling laws established in mammalian 
allometry arise from a natural similarity of animals and how they can be quite explicit on 
the evolution constraints through simple geometrical and clear physical conceptions. 

Acknowledgments 

The author is grateful to Per Christiansen and to Jayanth Banavar for their scientific 
interest to the current research and to the Ph.D. student Guilherme Garcia for helpful 
discussions of the biological literature. Financial support by the CNPq is also acknowledged. 



17 



Bone 


dimen sions 


reduc ed dim ensions 


bone mass 


muscle area 


LSR data 


di 


k 


k/di 


A. 


h*/di - 2 


2di + k 


hi 


buckhng 


bending 


humerus 


.3816 


.2996 


0.785 


0.763 


0.804 


1.063 


1.070 


0.927 


0.838 


radius 


.3868 


.2995 


0.774 


0.753 


0.802 


1.073 


1.084 


0.948 


0.850 


femur 


.3548 


.3014 


0.849 


0.843 


0.988 


1.011 


1.060 


0.816 


0.714 


tibia 


.3600 


.2571 


0.714 


0.764 


0.717 


0.977 


0.978 


0.926 


0.822 


Averaged 


.3708 


.2894 


0.781 


0.781 


0.828 


1.031 


1.048 


0.904 


0.806 


RMA data 


■di 


..li 


.. U/di 


A,- 


■ ..h/d-2 


2di + U 


hi 


(buckl) 


(bend) 


humerus 


.3860 


.3109 


0.805 


0.784 


0.806 


1.083 


1.083 


0.947 


0.862 


radius 


.4014 


.3210 


0.800 


0.787 


0.743 


1.124 


1.101 


0.959 


0.874 


femur 


.3599 


.3089 


0.858 


0.864 


0.976 


1.029 


1.071 


0.850 


0.799 


tibia 


.3654 


.2767 


0.757 


0.804 


0.731 


1.008 


0.998 


0.851 


0.781 


Averaged 


.3782 


.3044 


0.805 


0.810 


0.814 


1.061 


1.063 


0.901 


0.829 



Table 1. Testing of the mammalian long-bone similarity through the elastic-buckling and 
elastic-bending criteria. Experimental data by Christiansen on the mammalian dimension 
allometric exponents for i-bone diameter exponent di, length exponent li, reduced dimension 
exponent Aj, and bone mass exponent bi obtained by the least square regression {LSR) and 
the reduced major axis (RMA) regression methods. These are taken from Tables 2 in Refs. 
1^ , [p!5[] and |jl6[, respectively. The LSR data on b* are estimated here with the help of 



relation b* = Vibi , where rj (correlation coefficient) and 6, are corresponding data obtained by 
RMA regression. Predictions for the muscle-area critical exponents are given with the help 
of Eq.(16). Bone- averaged magnitudes are found as the mean values of the corresponding 
mammalian allometric exponents, e.g., d = Sj=i(ij/4. 
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Figure Captures 



Fig. 1. Evolution of the cross-section area for muscle fibers with body mass in the 
mammalian hindlimbs. Points: diamonds, circles, squares and crosses are experimental 



data taken from Fig. 3 in Ref. for, respectively, gastrocnemius, plantaris, deep digital 
flexors and common digital extensors. Arrows indicate the maximum muscle-area points 
achieved for a given mass. Solid line corresponds to regression Acm of these points along 
with their nearest neighbors, with Acm = 366 * M*^'^^. Dashes line is given for the isotropic 
scenario description with Aq = 29 * M'^/^. 



Fig. 2. Mammalian bone-dimension diagram: bone diameter against bone length. 
Points: A'79, B'83, B'92 and C'99 correspond to the overall-bone-averaged allometric data 
derived through the least square regression method by Alexander et at, Biewener, Bertran 
& Biewener, and Christiansen and reported, respectively, in Refs. P, [Tll , [l3| and 0]. Crosses 
correspond to the ESM p (rfg^"'^'^'-' = 3/8, l^^^''^^ = 1/4^ ao = bo = 1) and isometric scenario 
{do = lo = 1/3, ao = 2/3, 60 = 1) predictions; a-lines are due to the elastic-bucking and 
elastic-bending a-constraints given in, respectively, Eqs. (9) and (13) and estimated for the 
case of the critical-force exponent = 0.82 derived in Fig.l. The dashed areas indicate the 
elastic-bucking and elastic-bending criteria given, respectively, in Eqs. (10) and (14). These 
areas extend the corresponding a-lines by accounting of the 6-constraint equations within 
the experimental error for a(^P) = 0.82 ± 0.01 and 6("^p) = 1.05 ± 0.05 taken, respectively, 
from Fig.l and Table 1. 
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